Resilience of multi-photon entanglement under losses 
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We analyze the resilience under photon loss of the bi-partite entanglement present in multi-photon states 
produced by parametric down-conversion. The quantification of the entanglement is made possible by a sym- 
metry of the states that persists even under polarization-independent losses. We examine the approach of the 
states to the set of states with a positive partial transpose as losses increase, and calculate the relative entropy of 
entanglement. We find that some bi-partite distillable entanglement persists for arbitrarily high losses. 



I. INTRODUCTION 

Parametric down-conversion has been used in many ex- 
periments [1] to create polarization entangled photon pairs 
[2]. Recent experimental [3, 4] and theoretical [5-7] work 
has studied the creation of strong entanglement of large num- 
bers of photons. The states under consideration are entangled 
pairs of light pulses such that the polarization of each pulse 
is completely undetermined, but the polarizations of the two 
pulses are always anti-correlated. Such states are the polariza- 
tion equivalent of approximate singlet states of two potentially 
very large spins [8]. An application of the states for quantum 
key distribution has been suggested [5]. 

In any realistic experiment photons will be lost during prop- 
agation. It is therefore of great practical interest to analyze 
the resilience of the multi-photon entanglement under loss. 
A priori this seems like a very difficult task, because it re- 
quires the quantification of the entanglement present in mixed 
quantum states of high or actually even infinite dimensional- 
ity. However, the multi-photon states introduced in the above 
work exhibit very high symmetry - in the absence of losses 
they are spin singlets. The related symmetry under joint po- 
larization transformations on both pulses is preserved even in 
the presence of polarization-independent losses. This makes 
it possible to apply the concepts of 'entanglement under sym- 
metry' developed in Refs. [9-13] to the quantification of the 
multi-photon entanglement in the presence of losses. We cal- 
culate the degree of entanglement for the resulting states of 
high symmetry, as quantified in terms of the relative entropy 
of entanglement. We show that some (distillable) entangle- 
ment remains for arbitrarily high losses. 



II. SYMMETRY OF THE STATES IN THE PRESENCE OF 
LOSSES 

In the above-mentioned experiments and proposals a non- 
linear crystal is pumped with a strong laser pulse, and a three- 
wave mixing effect leads to the creation of photons along two 



directions a and b. To a good approximation the Hamiltonian 
in the interaction picture in a four-mode description is given 
by 



H 



b K (a{bt - alb{) + e- l K(a h b v - a v b h ). (1) 



The real coupling constant k is proportional to the amplitude 
of the pump field and to the relevant non-linear optical coeffi- 
cient of the crystal, and <j> denotes the phase of the pump field. 
Photons are created into the four modes with annihilation op- 
erators ah, a v , bh, b v , where h and v denote horizontal and 
vertical polarization. Note that both the modes and the asso- 
ciated annihilation operators will be denoted with the same 
symbol. In the absence of losses, this Hamiltonian leads to a 
state vector of the form [5, 6] 
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where r = nt is the effective interaction time and 
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In experiments the pump phase is typically unknown, and data 
is collected over time intervals much longer than the pump 
field coherence time. We will therefore consider the state p 
obtained from the state vector Eq. (2) by uniformly averaging 
over the pump phase <fi £ [0, 27r): 
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(n+1) tanh 2n r IV-X^-I- 
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The Hamiltonian H is invariant under any joint polarization 
transformation in the spatial modes a and b. That is, if one 
defines a = (ah,a v ) and b = (bh,b v ), then H is invariant 
under the joint application of the same unitary U from SU (2) 
to both vectors, a ^ U& and b i— > Uh. This invariance of 
H is inherited by the multi-photon states created through the 
action of H on the vacuum. This symmetry can be expressed 
as 
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for all U G SU(2), where V(U) = e mJ , and the real vector 
n is specified by U — e mCT / 2 , a denoting the vector of Pauli 
matrices. Here the angular momentum operator J can be writ- 
ten as J = J a + Jb- The components of J a associated with 
spatial mode a are given by the familiar quantum Stokes pa- 
rameters J a<x = (at_a+ — aLct-)/2, J a<y = (ajai — aj.a r )/2, 

and J Q:Z = (a\a h - oJo„)/2, with a± = (a h ± a v )/V2 
corresponding to light that is linearly polarized at ±45°, and 
a-i.r = [ah ± ia v )/ V2 to left and right-hand circularly polar- 
ized light. Analogous relations hold for spatial mode b. 

In the present work we are interested in the states created 
by H in the presence of losses. These losses will be modeled 
by four beam splitters of transmittivity r\ G [0, 1], one for each 
of the modes a,h,a v ,bh,b v , where the modes are mixed with 
vacuum modes. Explicitly, the operation corresponding to 
losses characterized by r] acting on a single mode a is given 
by 

OO 

£>) = X>>(Z«)t, (6) 

ri=0 

with L a n being given by 

t; = i=(i- ) ,)»oV a,B . (?) 

Vb! 

One can easily verify that these operators satisfy 

oo 

£(i-)ti« = l, (8) 

71 = 

required for trace preservation. In this paper we are interested 
in the situation where an equal amount of loss occurs in all 
four modes. We will denote the corresponding quantum oper- 
ation by 

£,=£° h ®£^®£j*<8>£|*. (9) 

It is not difficult to apply this loss channel to the state p of Eq. 
(4). However, the resulting expression is quite unwieldy, and 
quantifying the entanglement present in the state seems like a 
hopeless task at first sight. We will now discuss general prop- 
erties of the resulting state that allow a simple parametrization 
and as a consequence the determination of its entanglement. 

In the absence of losses, all components of the state created 
by the action of H have an equal number of photons in the a 
modes and in the b modes, since photons are created in pairs. 
The state vector \ip) of Eq. (2) is a superposition of terms cor- 
responding to different total photon numbers. For any given 
term we will denote the number of photons in the a modes by 
a = an + a v , where ah is the number of photons in mode ah 
etc. Analogously, the number of photons in the b modes is de- 
noted by (3 = (3h + (3 V . The relative phase between terms with 
different values of a or (3 depends on the pump phase <f>. The 
corresponding coherences in the density matrix are removed 
when averaging over the pump phase. 

Losses lead to the appearance of terms with a ^ (3. The 
state p 1 = C v (p) after losses now has the form 

oo 

p'= ]T P(a,f3)p^\ (10) 

a, (3=0 



where P(a, (3) is the probability to have photon numbers a 
and (3 in the a and b modes respectively, and p( Q ^) is the 
corresponding state. In the state before losses, the terms 
p( a , a ) are maximally entangled states (for a ^ 0), denoted 
by IV'-XV'-I in the notation of Eq. (3). Losses reduce this 
entanglement, but do not make the state become separable, as 
will be seen below. 

The state vector \ah, a v )\/3h, Pv) corresponds to a spin state 
vector \j a , m a )\j b , m b ) with j a = (a h + a v )/2, m a = (a h ~ 
a v )/2,j b = {I3 h + f3 v )/2,m b = {fi h - (3 v )/2. Note that 
in this representation a single photon corresponds to a spin- 
1/2 system. A state with fixed photon numbers a and j3 thus 
corresponds to a state of two fixed general spins j a = a/2 
and j b = /3/2. 

The key feature of the lossy channel C„ of Eq. (9) is that it 
does not destroy the symmetry described by Eq. (5). We have 
that 

V{U)C v {p)V(Uy =C n (p) (11) 

for all losses rj and all U G SU(2). To sketch the argument 
why this symmetry is retained we will resort to the Heisenberg 
picture. Polarisation-independent loss in the a modes can be 
described by the map 

ana' = y/rja + y/l - r/c, (12) 

where c = (c/j, c v ) is a vector of unpopulated modes that are 
coupled into the system due to the loss. Applying U G SU (2) 
to a' gives 

a" = Ua' = VWa + V 1 - r lUc. (13) 

On the other hand, applying first U and then the loss operation 
gives 

a" = 0)f/a+\/r^c, (14) 

in which the last term is different. However, this term just cor- 
responds to a coupling in of unpopulated modes with a coef- 
ficient \/l — rj. The resulting lossy channel is invariant under 
the map c i— > Uc, since these modes are unpopulated. This 
implies that the state after application of the loss operation C v 
has the same symmetry as before. Note that for this argument 
to hold, the amount of loss in the a and b modes does not have 
to be the same, since the transformations are applied indepen- 
dently to each of a and b. However, within each spatial mode, 
losses must be polarisation insensitive. 

The identification of the above symmetry dramatically sim- 
plifies the description of the resulting states. The most gen- 
eral state p( a '' 3 ) with fixed value of a and (3 for which 
V(U)p( a >&V(Uy — p( a >@^ for all U G SU (2) is of the form 

P^= Xf ^fif'fl (15) 

j = \ja-jb\ 

where j a = a/2,j b = (3/2 [13], essentially as a consequence 
of Schur's lemma [14]. Here, the p[ a '^ form a probability 



3 



distribution for all (a, (3) in the allowed values for j. In turn, 
flj is up to normalization to unit trace a projection onto 
the space of total spin j (for fixed j a = a/2,jb = (3/2). That 

is, = l^' m /(2j + 1), where tf' P) is equal to the 

identity when acting on the space labeled by a, (3, and j, and 
zero otherwise [13, 15]. 

As an example, let us consider the case with exactly one 
photon in each spatial mode, i.e., a = (3 = 1. Then 
there are just two terms in the expansion of Eq. (15), pro- 
portional to fig and Qj 1 ' 1 . The state Slg 1,1 -* is the projector 
onto the two-photon singlet state with state vector ((a^J, — 
albjj j \/2)|0) , while fi^ 1 ' 1 ' is the normalized projector onto 
the spin-1 triplet. The trace condition (Iq + Mi = 1 
means that the set of all invariant states p^ 1 ' 1 ' is characterized 
by just one parameter. Note that the most general state with 
exactly one photon in each spatial mode would be character- 
ized by 15 parameters. 



III. QUANTIFYING THE ENTANGLEMENT 

In order to quantify the entanglement in a given physical 
situation, one has to determine the coefficients P(a, (3) of 

Eq. (10) and Pj of Eq. (15), which may be calculated 
from the polarization dependent photon counting probabilities 
p(ah, a v , f3h, f3 v )- These in turn can be determined by explic- 
itly applying the loss channel of Eq. (9) to the state p of 
Eq. (4). One finds 



rameter Pq~ is given by 



p(a h ,a v , (3 h , (3 V ) 

oo 

X E 



(cosh(Kt)) 4 a h \a v \l3 h \p v \ 
((l-r;) tanh(^)) 2 ( m + n )( 7 n!) 2 (n!) 2 



(16) 



(m - oth)\{m - (3 v )\(n - a v )\(n - (3 h )V 



where itiq = max(a/,,/3„) and no = max(a„,/3/,). The prob- 
abilities P(a, (3) are obtained by summing this expression 
over all ah,a v ,f3h, (3 V with at + a v = a and (3h + f3 v = 13. 

The coefficients p™ may be written as linear combi- 
nations of the p(ah, a v , f3h, (3 V ) via the Clebsch-Gordan co- 
efficients [14] by means of the standard procedure of 'cou- 
pling spins'. Polarization-sensitive photon counting in the 
spatial modes a and b corresponds to the basis spanned by 

the \ j a , m a ) |jb, rrib), while the p™ and Q^'^ are defined 
in terms of the total spin, corresponding to the label j. Since 



the p" 1 ' characterize the normalized state p^ a, ^\ they only 
depend on the relative probabilities of the different values of 
a/j, a v , (3h, (3 V for given a and (3. Eq. (16) then implies that 
they depend on the interaction time t and the transmission r\ 
only via the combination £ = (1 — if) tanh(KTj) <E [0,1], which 
ranges from zero for perfect transmission (or, less interest- 
ingly, zero interaction time) to one in a limit of complete loss 
and infinite interaction time. 

For example, for a = (3 = 1, the single independent pa- 



(i>i) 
Mo 



1 - -( P (l > 0,l,0)+j)(0,l > 0,l))/P(l,l) 



(17) 



where as before P(l,l) = p(l, 0, 1, 0) + p(l, 0, 0, 1) 
p(0, 1,1,0)+ p(0, 1, 0, 1). This gives 



Mo 



(i + e 2 /2)/(i + 2f). 



(18) 



To quantify the entanglement present in the total state, one 
can proceed by considering each pt"'^ separately. There is 
no unique measure of entanglement for mixed states. Instead, 
there are several inequivalent ones, each of which is associ- 
ated with a different physical operational interpretation [16]. 
The relative entropy of entanglement [18], which will be em- 
ployed in the present paper specifies to which extent a given 
state can be operationally distinguished from the closest state 
that is regarded as being disentangled. The relative entropy of 
entanglement of a state p is defined as 



Er(p) 



inf Seller), 



(19) 



where 5(jo||cr) = tr[plogp — ologer] denotes the quantum 
relative entropy of the state p relative to the state a. Here V 
is taken to be the set of states with positive partial transpose 
[19] (PPT states). This set of states includes the set of separa- 
ble states, but in general also contains bound entangled states 
[20] . The relative entropy of entanglement is an upper bound 
to the distillable entanglement [16], providing a measure of 
the entanglement available as a resource for quantum infor- 
mation purposes [21]. 

The symmetry of the states dramatically simplifies the cal- 
culation of the relative entropy of entanglement. As follows 
immediately from the convexity of the relative entropy and 
the invariance under joint unitary operations, the closest PPT 
state can always be taken to be a state of the same symmetry 
[10, 13]. Hence, the closest PPT state is characterized by the 
same small number of parameters. For simplicity of notation, 
we will denote the subset of state space corresponding to spe- 
cific numbers a, (3 of photons as (a, /3)-photon space. In the 
(1, 1) -photon space let us denote the closest PPT state as 



(20) 



Forming the partial transpose of this state, and demanding 
that the resulting operator be non-negative, gives the condi- 
tion Co 1 ^ < 1/2- In this simplest space, all symmetric states 
lie on the straight line segment p$ £ [0, 1] with the PPT 
region extending from the origin to the midpoint (see Fig. 1). 

In general, for higher photon numbers a and (3, the 
set of symmetric states are represented by a simplex in a 
(min(a, (3) + 1) -dimensional space, the coordinates of which 

are denoted by (ij . In turn, the PPT criterion gives rise to a 
number of linear inequalities, such that the set of invariant op- 
erators with a positive partial transpose corresponds again to 
a simplex. The intersection of the two simplices corresponds 



4 



(1/2,1/2) 
— -x — 




(0,1,0) £ ^ (1,0,0) 

Or. 




FIG. 1: The simplices of symmetric states for the cases of (a, j3), 
a = j3 — 1,2,3, respectively. The equilateral triangle has been 
marked with contour lines, on each of which one of the parameters 
is constant. The set of PPT states is indicated by the grey line seg- 
ment in the top graph, the shaded area of the (2, 2) triangle and the 
filled polygon which obscures the fi^' 3 ^ = 1 vertex in the (3, 3) 
tetrahedral space. In all graphs only the projector of highest spin is 
within the PPT set. For all three cases, the set of all possible down- 
conversion states is a curve ending at the boundary of the PPT set, 
shown by the solid black line for the (1,1) space, and by the dotted 
curves for the (2,2) and (3,3) spaces. The position of the state on the 
curve is determined by the parameter £ = (1 — r)) tanh(ftt). 



to the invariant PPT states, and the coordinates are denoted by 
C C«,/S) [22] _ 



The situation with a = (3 = 1,2, 3 is depicted explicitly in 
Fig. 1 . The simplex corresponding to symmetric states, char- 
acterized by the condition that the form a probability 
distribution, is in these three cases a straight line segment, 
an equilateral triangle, and a regular tetrahedron respectively. 
The vertices of the simplex represent the normalized projec- 



tors QS 



,0) 



States in the interior of the simplex are convex 



combinations of all the allowed projectors. The PPT set with 
the same symmetry is clearly marked. 

Fig. 1 also shows the curves traced by the down-conversion 
states when they are subject to loss. As discussed above, the 
position of the states on the curve is determined by the single 
parameter £. For perfect transmission corresponding to r] = 1 
the quantum state in an a = (3 photon space has [j,q*'^ = 1 
for all values of t, corresponding to maximal entanglement. 
As losses are increased the state migrates towards the PPT 




1.0 Tl 



FIG. 2: Lower bounds to the relative entropy of entanglement for 
down-conversion states with initial average photon numbers of 0.5 
(solid), 1 (dashed), and 3 (dotted line) subject to loss, evaluating the 
sum of Eq. (21) up to a truncation of a, /3 < 5. This gives a good 
approximation to the total entanglement for average photon numbers 
before loss up to about 3. 



boundary. It is an important immediate consequence of Eq. 
(18) that for all losses rj > 0, the number /4 ) Q ' Q ' ) is always 
greater than 1/2 for any finite t and for all a. For any fi- 
nite t, -^l/2as£^l (which corresponds to a limit 
of zero transmission time and infinite interaction time). This 
holds true for (a, a) = (1, 1), but also for higher values of a: 
the state remains outside the PPT set for any non-vanishing 
t and for arbitrarily high losses. Therefore, the above results 
show that there is always some entanglement in the down- 
conversion state, as quantified in terms of the relative entropy 
of entanglement. As a corollary, which one can already infer 
from the lowest dimensional subspace, (a, a) = (1, 1), there 
is actually distillable entanglement in the down-conversion 
state, regardless of how lossy the transmission from the source 
to the detector. 

We now proceed to quantify the entanglement in the states 
more explicitly. Since Er is convex and the set of symmetric 
PPT states is convex, finding the closest state a amounts to 
solving a convex optimisation problem. For different values 
of a, (3 the quantities S{p^ a '^'\\(j^ a '^) have been evaluated, 
where cr( Q ' 3 ) denotes the PPT state which is the unique global 
minimum in the convex optimization problem, i.e., the PPT 
state closest to the down-conversion state. For generic states, 
this optimization problem would still be convex, yet, the di- 
mensionality of state space grows as (a+l) 2 (f3+l) 2 — 1. The 
symmetry dramatically reduces the dimensionality of the con- 
straint set to searched to min(a, (3), and thus makes the quan- 
tification of the entanglement a feasible task. For instance, 
for a state with three photons on each side, one has to con- 
sider only three objective variables instead of 255. The total 
relative entropy of entanglement is given by the expression 



Er{ P )= £ P{ a j3)E R (p^). 



(21) 



a,/3=0 



The average photon number before loss N is related to the 
interaction time t as N = 2sinh 2 (Ki). The average photon 
number after loss is n = i]N. Fig. 2 shows the relative entropy 
of entanglement calculated as described above for N = 0.5, 
N = 1 and = 3. One sees that significant entanglement 
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remains even for substantial losses. 



IV. CONCLUSIONS 

We have shown how symmetry considerations make pos- 
sible the quantification of entanglement for states produced 
by parametric down-conversion and subject to losses. The 
resilience of the entanglement of these multi-photon states 
under photon loss makes them an excellent system for the 
experimental demonstration of entanglement of large photon 



numbers [4] and good candidates for quantum communication 
schemes [5]. 
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